Gapped Heisenberg spin chains in a field 
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We consider the fully anisotropic Heisenberg spin-1/2 antiferromagnet in a uniform magnetic 
field, whose ground-state is characterized by broken spin rotation symmetry and gapped spinon 
excitations. We expand on a recent mean-field approach to the problem by incorporating fluctuations 
in a loop expansion. Quantitative results for the magnetization, excitation gap and specific heat are 
obtained. We compare our predictions with new DMRG and exact diagonalization data and, for 
zero field, with the exact solution of the XYZ spin chain from the Bethe Ansatz. 

PACS numbers: 75.10.Jm 



I. INTRODUCTION 

For many years now, low-dimensional quantum sys- 
tems have been a hotbed of theoretical and experimental 
novelties. Quasi one-dimensional spin systems have by 
now become fertile laboratories for the investigation of 
many interesting phenomena stemming from strong cor- 
relations, an important example being charge fractional- 
ization of physical excitations. Overall, spin-1/2 Heisen- 
berg chains count among the most important paradigms. 
This stems first from the existence of many exact results 
on their ground-state, excitations and thermodynamics 
in a variety of cases, but also on the growing number 
of experimental situations in which the effective theory 
boils down to a one-dimensional exchange Hamiltonian. 

Exact results on quantum spin chains (by which we 
mean solvability at the lattice level) are constrained to 
either easy-plane anisotropy in nonzero longitudinal field 
(so spin projection along the field axis is a good quantum 
number, and spinon excitations are well-defined) (for a 
full exposition of this well-known subject, see e.g. 0l3l)i 
or zero field and general anisotropy^. Recently, how- 
ever, there has been renewed interest in the study of spin 
chains in t he m ore general case of transverse field and 
anisotropy |4l5l lal7^ , where all symmetries are broken and 
exact solutions are out of reach. Most notably, the main 
motivation for the present study originates from exten- 
sive work done on the field-dependent ground-state prop- 
erties and excitation dynamics of quasi-one-dimensional 
anisotropic antiferromagnets like Cs2CoCl4^. 

We will consider the Hamiltonian of the generally 
anisotropic Heisenberg chain in a field given by 

N 
3=1 

where the spin-1/2 operators S*" obey the commutation 
relations Sf,S^, — iSjj'e°''^'^Sj . This model shares 
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many features with somewhat simpler and better un- 
derstood exchange Hamiltonians like the XY— or XXZi*^ 
models. As far as the XY model is concerned, many ex- 
act results are known in the presence of a field^*^, like 
the existence of an Ising-class disordering transition at 
a critical value of the field. General space- and time- 
dependent correlation functions can be obtained by map- 
ping to gapped noninteracting fermion3ifi*ii*iSiiii4. 

On the other hand, it is well-known that the isotropic 
antiferromagnetic Heisenberg chain, and more generally 
the XXZ {Jx = Jy y^ Jz) chain in a field parallel 
to the anisotropy (with —1 < Jz/Jx < 1), is a gap- 
less system with no long-range order and an unbroken 
t/(l) symmetry. Once again, many results are known on 
the thermodynamicai^ while recent progress has been 
achieved on correlation functions both on the latticei^ 
and in the scaling limilsi^. 

For the fully anisotropic model in a field iQJ, how- 
ever, precious little is known. In zero field, the model 
becomes exactly solvable'^, meaning that thermodynamic 
properties are computable^. Correlation functions in zero 
field are still unknown in general (with the exception of 
the Baxter-Kelland formulaiSiiS for the one-point func- 
tion giving the staggered magnetization), although some 
special cases have been studied^SiSi. The distinguish- 
ing feature of Q is that one can find two independent 
origins for a gap: first, in zero field, the very presence of 
structural anisotropy breaks the symmetry right from the 
start. One could also break the symmetry by applying a 
field transverse to the anisotropy (like in the case of an 
XXZ magnet with field along the y direction). In both 
of these cases, long-range Neel order develops. However, 
the noncommuting nature of the field with respect to the 
anisotropy and the fact that the field strength is tunable 
lead to the existence of quantum critical behaviour in 
the system: for large enough fields, fluctuations become 
strong enough to disorder the ground-state, leading to an 
Ising-class critical point at a field H = He (whose value 
depends on the anisotropy) . This holds true for a general 
anisotropic chain, studied in )17l |. in which bounds were 
obtained for the critical field, and rigorous results were 
derived on classical ground states for particular values 
of the field (meaning that the ground-state becomes a 



two-sublattice Neel-type state, with nonzero long-range 
antifcrromagnetic order). 



Exact results being inaccessible for now, one must rely 
on alternative methods to compute quantities of interest 
like the magnetization along the field, transverse order 
and objects like the dynamical structure factor. Close 
to a gapless, solvable point, one can rely on approaches 
like bosonization to obtain various scaling results {ej/. 
weak field and/or small anisotropy limits studied in|a). 
However, the presence of a gap justifies the hope that 
a mean-field scheme might work out relatively well as a 
quantitative theory provided the correct starting point is 
identified. This was initiated in [J| for the XXZ chain 
in a transverse field, where the phase diagram and vari- 
ous scaling results were obtained. Correlation functions 
were afterwards studied in Q using the equivalence to 
the XY model in a field, and through the use of per- 
turbative RG shemes based on sine-Gordon form factors. 
Some of these results are in fact immediately general- 
izable to the XYZ chain in a field that we purport to 
study. One thing that was missing, however, was an es- 
timate of the strength of fluctuations, which is of partic- 
ular importance in view of potential quantitative com- 
parisons with experiments. For low fields, the mean- field 
approach failed by construction to properly describe the 
restoration of the U{1) symmetry of the XXZ chain. The 
mean-field staggered magnetization, for example, failed 
to vanish at zero field, in clear contradiction with the 
exact result. The field dependence of the gap also failed 
for small fields, disagreeing with both experiments and 
numerical results using the Density Matrix Renormaliza- 
tion Group (DMRG). One might ask whether DMRG or 
other quasi-exact methods would suffice to describe all 
aspects of experiments, but its main drawback is that 
it is only very useful as a tool for investigating prop- 
erties of the ground-state and lowest-lying states, and 
is not yet an economical tool for dynamical quantities 
like spin correlation functions, measurable using neutron 
scattering. Our general idea is therefore quite simple: to 
carefully compute fluctuations around mean-field to two 
loops, providing an improved setup from which values for 
the thermodynamic (and eventually dynamic) properties 
of the system, and quantitative estimates of the precision 
of the mean-field theory can be obtained through the use 
of intensive numerics on the resulting equations. 



This paper is organized as follows. In section If, the 
basic mean-field setup is presented. Fluctuations around 
mean-field are computed in section III, and used to de- 
termine corrections to physical parameters, the specific 
heat, and the excitation gap. Section IV presents a com- 
parison of our approximations in zero field with the exact 
results available from the Bethe Ansatz, and conclusions 
are collected in Section V. 



II. SETUP 

Written as it is, (Q) is not tractable in any simple way. 
To make progress, we move from a spin basis to a sim- 
pler fermion one by using the Jordan- Wigner transforma- 
tion 5'! = cjcj - 1/2, S*] = e*'^Sill4='=c], under which 
the Hamiltonian gets transformed to a theory of inter- 
acting fermions in the presence of a superconducting-like 
Cooper pairing term: 

N J j_ 

+Uc]c, - l/2)(c]+ic,+i - 1/2) - H,{c]c, - 1/2) (2) 

where J± = [J^ ± Jy)/'2.. We impose periodic bound- 
ary conditions along the lattice on the spin operators, 
so iS^^jY = S*". The fermion operators are thus taken 
to have antiperiodic boundary conditions. We shall re- 
strict the exchange parameters to the domain Jx^ Jz^ 
0,Jx > \Jy\ > 0, therefore concentrating on the antifcr- 
romagnetic regime. Note that (Q is not quite the most 
general anisotropic exchange Hamiltonian in a field that 
we could have started from: pointing the (field, or equiv- 
alently the ) anisotropy in a completely general direction 
would entail terms of the form S°'S^^^ with a j^ (3. The 
absence of cc, z and y, z exchange terms is here required 
to ensure a local fermionic theory. 

The interacting fermion Hamiltonian (0) is not ex- 
actly solvable because of the simultaneous presence of 
the pairing-like term J_, originating from the transver- 
sality of the magnetic field and x, y anisotropy, and the 
external magnetic field. We therefore use the approxima- 
tion scheme proposed in j^] and elaborated in ,5] . Let us 
recall here the basics of this setup. First, the fermionic 
interaction term is decomposed via the (site independent) 
mean-field parameters 

■^0^ {{cjCj -l/2))o, /Co = (c]+iCj)o, ■^o = (cj+iCj)o(3) 

corresponding to the magnetization, kinetic and pairing 
terms respectively. The angular brackets with subscript 
signify thermal averages using the mean-field states. The 
mean-field Hamiltonian becomes an exactly solvable free 
fermion theory (corresponding to an effective anisotropic 
XY model in a field), 

^ t A ■ 



J=l 



-h{c]cj - 1/2) + eo 



.,., .,_,,.. (4) 

in which eo = Jz{-Ml + /Cq - Vq), and 

f-J+-2J,/Co, A^J^+2J,Vo, h = H,-2J,Mo.i5) 

The parameters t, A and h are thus simply in- 
terpreted physically as effective bandwidth, effective 



superconductivity-like pairing, and effective chemical po- 
tential. The ground-state is therefore BCS-like and con- 
tains all particle numbers. Performing a Fourier trans- 
form to momentum space using Cj = — ^ 12n=i 6''^""'cfc„, 
kn = ^{n — 1/2), Cj+N — —Cj and choosing TV even for 
convenience, we can write 



N/2 , 

■^-^ \ I t cos kn — h j A sin fc„ 



Cfc„ + eo (6) 



using the Nambu spinor cJ, — { cl c^k„ ] ■ The theory 
can easily be rediagonalized as 



N/2 



Hmf = ^c'l^a''u{kn)c',. 



k„ +eo 



(7) 



with Bogoliubov-de Gennes transformation parameters 



Cfc — BfcCfc, Bfc 



a+(fc) ia-{k) 
ia^{k) Q;+(fc) 



a±{k) = 



^/2 



1± 



t cos k — h 



-,1/2 



;(fc) 



and a±{~k) = ±a±(fc). The dispersion relation is 



u;{k) = [{t COS k-hf + A^ sin^ fc] 



1/2 



(8) 



(9) 



(10) 



and makes manifest the presence of the gap in the the- 
ory. Now, in terms of the new fermions, we sit always 
precisely at half-filling, and the ground-state is charac- 
terized by a fully filled lower band. Particle-hole excita- 
tions correspond to spinon-like excitations in the original 
spin model: exactly as in the XY model, due to energy- 
momentum conservation two-particle processes dominate 
below the critical field, whereas above it one-particle pro- 
cesses give the strongest contribution to correlation func- 
tions. 

The thermodynamic properties of the model at the 
mean-field level are easily obtained from the mean-field 
free energy per site at inverse temperature j3, 



/o 



-1 



InZ 



M_F, 



Explicitly, in the limit N 



Zmf = TTe-'^"""'. 



cxD, we get 



/o 



eo 



-2 r dfc , , (3uj(k) 
—— I T^mcosh — - — . 



P 



27r 



(11) 



(12) 



The mean-field setup is completed by requiring that the 
following self-consistency equations be satisfied: 



'dk 



Mo + l/2 = y"g [l-{l-2np{u(kmal(k)-a'_{k))\ , 

/Co- / — cosA: [l~{1^2np{uj{k))){a\{k)^al{k))\ , 
Jo ^''' 

(13) 



Hdk 
Po = - / — sinfc(l-2n0(w(fc)))a+(fc)a_(A:) 
Jo 27r 



with nfj{a) 



(l- 



cl3a\ 



~^ . These are easily solved numer- 



ically with a simple iterative procedure. One important 
consequence is that all mean-field parameters, and there- 
fore all effective potentials appearing in the mean-field 
Hamiltonian, acquire a temperature dependence. The 
average energy per site in the mean-field approximation 
is itself simply given by 



Mfc ,,, , (3uj{k) 
eo= eo- / ^w(k) tanh — - — 

tK-Q + ZSPa-hMa + eo. 



(14) 



The mean-field specific heat is defined for fixed external 
field iJ^, c-^^l 



dT 



H, 



III. FLUCTUATIONS AROUND MEAN-FIELD 

The fluctuations around mean-field can be described 
using standard pcrturbative expansions by reinstating 
the original fermionic interaction term in the Hamilto- 
nian Q . In detail, the perturbation that we have to deal 
with reads 



H^ 



J: 



FL 



N 
ni,n2,"3 = l 



^"3 



H+"3 



Cfc„ 



efe„ :(15) 



in which the symbol :: signifies that self-contractions are 
absent in perturbation theory, and where we introduce 
the notation fc° = ^n. We are being pedantic with 
the momenta, since we keep all finite N corrections in 
the computations (this turns out to help the numerics 
perceptibly, since we can realistically plot our functions 
for N only of the order of a few hundred). 

Since we are interested in equilibrium thermodynamic 
quantities, the pcrturbative computations are most eas- 
ily tackled using the imaginary time (r = it) formal- 
ism. First, we define the Matsubara Green's functions 
as appropriate averages over the mean-field theory (rep- 
resented everywhere by (...)o). We need two types of 
functions, one representing normal correlations and the 
other, pairing correlations: 

Go{kn,-iT)6k^^k^, = -(Tt- |cfe„(-ir)4^, (0)|)o 
J"o(fc«,-«T)4„+/c„,,o = (Tr {cfc„(-ir)cfc„, (0)})o. (16) 

The most convenient representation, exhibiting the poles 
in the normal and Cooper channels, is obtained after 
Fourier transforming to frequency space (we use the def- 

initions fiiujc) = jyTe'^-^fi-ir), w„ = f{a + 1/2), 
ag Z) 



go{k,iuJa) = - 



_{k) 



-(k) 



iuja — uj(k) iLUa+l-u{k)' 



J-o{k,iuJa) — 2a+{k)a-{k) 



iuj{k) 



^2 -f w2(fc)' 



(17) 



All physical quantities of interest can be computed using 
these functions, in particular the unperturbed mean-field 
parameters: 



1 ^ 



[3N 



n=l q£ 

N 



ICq = — y^ y^ cosknGojkn, iujg), 

n=l agZ 



N 



Vo 



(3N 



^^smknJ^o{kn,iuJa). (18) 



in which all correlators are now evaluated at mean-field 
level. To first order in perturbation theory, fluctuations 
do not correct the mean-field Green's functions (by def- 
inition). The first corrections appear at second order, 
where for each type of Green's function 12 different dia- 
grams can be written down. Summing these after chop- 
ping external legs provides a computation of the self- 
energies to two loops. After a relatively lengthy compu- 
tation, the self-energies to second order (as identified by 
the superscript) can be written down as (after summa- 
tion over internal Matsubara frequencies, only the sum 
over the two internal momenta ki , /c2 remains) 



The full Green's functions in imaginary time are de- 
fined as 

g(fc,-*r) = -(T,{cfe(-jr)4(0)}), 

nk, -it) = (T, {cfc(-zT)cfc(0)}) (19) 

where the angular brackets denote a thermal average over 
the full Hamiltonian H = Hmf-'tHpl. The perturbation 
expansion can then be formally written down as the series 
(here in the normal Green's function) 

°° f— 11" f^ 

g{k, -it) = - V ^-^ / dn...dTn X 

x(T, {ckHT)HFLi-^Tl).■■HFLi~^Tn)cliO)})o (20) 



4j: 



Y!-^\k,iuj) = -| ^ TAr(fc,Zw|fci,fc2), 



AJ: 



T.f{k,iuj) = ^^Y.Ts{k,^uJ\k,,k2), (21) 



fci,fc2 



where the "normal" and "Cooper" self-energies both re- 
ceive contributions from normal and Cooper channels. 
Explicitly, we write those contributions as 



T^(fc,zc.|fci,fe)- Y. {^~5)[C^^'',^,^{k\k,,k^)h{k\k^,k2) + C^,^,'l,j2{k\k,,k2^^ 



Ts{k,iuj\ki,k2) 



(Ti (72 Cr3 
CTl (72 (73 



Ni 



(1-5) [CZl.,mi,k2)fi{k\ki,k2) + C^^f,3/2(fc|A:i,fc2)] D,,,,,,{k,tLu\ki,k2) + STI, {22) 



and the external frequency itj-dependent part is 

[n;3(o-3w(fc„J) - n/3(cr2w(fc„j+„J][n^(o-2w(fc„i+„2) - cr3w(/c„J) -|- n^(o-iw(fc„j+„))] 



^(Ti(T2(T3 \^n 5 ^^|"'ni 7 "-r? 



iuj - aiU!{km+n) + <J2U}{km+n2) - o-3w(fc„2) 

I 



(23) 



where n^{a) = (e''° - 1)"^ The terms T^g are finite N 
corrections which are detailed below. These have been 
explicitly kept in the numerics, where they noticeably 
accelerate convergence towards the thermodynamic limit. 
The parameter 



rameters, namely 



S - (5<T20-3(5ij(fc„j+„2),c^(fc„2) 



(24) 



comes from the double-pole contributions in the inter- 
nal Matsubara frequency summation. It can of course be 
put to zero in the strict N ^ oo limit. The momentum- 
dependent parameters in (|22|) are given by various combi- 
nations of the Bogoliubov-de Gennes transformation pa- 



C'aiJ^as = crio-2a+-(fc„i+„)a+_(fc„,+„Ja^,,(fc„J, 

Cfi'i^o-3 = (Ji(J2(J3a+-{kn,+n)a+-{kn,+n2)a+-{kn2), 
(^a,a2% =Cr3al,{kn,+n)al^{kni+n2)a+-{kn2) (25) 

where a-| (fc) — a+{k)a^{k), and we have suppressed 

the arguments (k\ki,k2) for clarity. The momentum 
weight factors are trigonometric functions given by (we 
remind the reader that k^ = •^n; moreover, since we 
sum over the full Brillouin zone, the Umklapp terms are 



explicitly included in our expressions) 

f (I. \i. I. \- (cosfcO^ -cosfc,",^_„)^ 

J2['^n\kni J n^n2 j ~ (COS fcn^-n ~ '^OS fc„^ ) X 

x(cosfc°^_„-cosfc°^^^^+i). (26) 



expressions (|22() are finally given by 






The finite N correction terms included in the self-energy in which (for A,B,C^ N, S) 
I 



(27) 



V 



ABC 



ftikl+k) 



/^(fcl + fc) 

iu> — uj{ki + k) ' iuj + uj{ki + k) 



(/f(fc2)/t(fc2)+/^(fc2)/^(fc2))nM^(fc2)(l-nMc.(fc2))) (28) 



and where we have used the notations 

f^ik)^alik), /|(fc) = T*a+-(fc). (29) 

It is evident that the self-energies manifestly obey the 
same symmetries as the Green's functions, namely 



self-energy matrix 



J^jS[{k,iuj) — J^]\j{—k,iuj), 
T,s{k,iuj) = -T,s{-k,iuj), 
Ss(fc,-iw) = ^s{k,iuj)- 



(30) 



At this stage, a perturbative (partial) resummation can 
naturally and easily be performed using Dyson's equa- 
tion. We start by defining the basic propagator matrix 

go{kM=( %^)C'"\ ?S''"\V (31) 

^ ' \-To(k,tuj) -yQ(k,-iuj) J ^ ' 

In the same way, we can arrange the self-energies into the 



J 



S(fc,iw) 



YiN{k,iuj) I]s(fc,iw) 
— Es(fc,ia;) —T,pf{k,~iLL)) 



(32) 



The basic resummation is then performed by noting that 
the full propagator matrix fulfils the Dyson equation 

Q{k, iuj) = Qo{k, iuj) + Qo{k, iuj)'E{k, iuj)Q(k, iuj), (33) 

in turn leading to the full matrix propagator 

g{k, iu;) = [g^^k, iuj) - S(fc, iuj)] "^ . (34) 

In terms of mean- field Green's functions and self-energies, 
this reads 



g{k,iuj) ^ D-^ 



go{k,iuj) + Dct{QQ{k,iuj))T,i^{k,~iu;) Ta[k,iuj) + Dct{Qfi{k,iuj))T,s{k,iuo) 
— JFo(fc, ibj) — Det(t;o(fc, iu}j)'Ss{k, iuj) —Go{k, —iuo) - Det(tJo(fc, iuj))T,]y{k, iu) 



(35) 



r 



with denominator 



A. Corrections to mean-field parameters 



£'=1 -Tr(ao(fc, iw)S(fc, iw))-f Det(go(fc, iw)S(fc, iuo)) 



(36) 



As a first application of our formulas, the modifications 
to the mean-field parameters can be computed using our 
full Green's functions. We define the general (off mean- 
field) averages 



where the operators Tr and Det operate in the Nambu 
space, and where Det(Po(fc,i^)) — —[u'^ +Lo'^{k)\^^ . The 
two-loop approximation then corresponds to replacing 
the exact S(fc, iuj) by the two-loop result Yp''{k^ iuj) (|21|l . 



M = ((c]c,-l/2)), IC = {c]^,c,), r^{c,+^c,), (37) 



Magnetization vs. field: DMRG and MF+fiuctuations 



Kinetic term vs. field: MF and MF+fiuc 




mean field 

fiuctuations 
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DMRG N=40 n 

DMRG N=50 ■ 

DMRG N=60 o 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 



FIG. 1: Magnetization as a function of field, for different 
anisotropy cases, computed using the basic mean-field setup, 
and with fluctuations. The data points are the DMRG results. 



FIG. 2: Kinetic term as a function of field, for different 
anisotropy cases, computed using the basic mean-field setup, 
and with fluctuations. 



which can (in complete analogy with the pure mean-field 
ones, see (I18fl ) be obtained via the summation rules 



Pairing vs. fieid: MF and MF+fiuc 



1 ^ 

1 ^ 
^ = l^M^ y^COSfcng(fcn,ZLJa), 

N 

'^ = J^J2Y1 sinA:„^(fc„,ic^o). (38) 

^ n=l aSZ 

Plots of these quantities as a function of field for various 
values of the anisotropy parameters are given in figures 
(1-3). In the case of the magnetization, it can be seen 
that the inclusion of fluctuations within the framework 
has vastly improved the predictions of the theory all over 
the parameter space. For example, the relative error at 
H, = 1 for iJa:,Jy,Jz) = (1.0,0.25,1.0) has gone from 
12.5 percent for the mean-field versus DMRG result, to 
1.9 percent for fiuctuations versus DMRG. 



B. Energy per site and specific heat 

The mean energy per site can be simply obtained per- 
turbatively either directly or through a linked cluster-like 
computation as 

e = ea+t5IC + A5V-hSM + 



1 



4:NP 



^ Tr[g(fc,iw„)S(fc,iw„)] (39) 



in which SIC = /C — /Cqi etc. Substituting the two- loop 
results we have derived for the matrix Green's func- 
tion and self-energy, we can compute the mean energy 
and more interestingly the specific heat as a function 
of field and temperature. In figure (01, the XXZ-like 
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FIG. 3: Pairing term as a function of field, for different 
anisotropy cases, computed using the basic mean-field setup, 
and with fluctuations. 



case with anisotropics (1.0, 0.25, 1.0) is considered. Since 
this system is gapless at zero field, this is the "worst- 
case" scenario considering that our mean-field theory 
starts with a nonzero gap even at zero field. Both the 
straight mean-field, and the fiuctuations-corrected val- 
ues are plotted. The gap computed using fluctuations is 
lower than the mean-fleld one, as can be seen from the 
low-temperature onset of the specific heat. For external 
field Hz = 1.8, above the transition, the gaps are essen- 
tially the same, and so are the specific heat curves. In 
figure (O, the differences between the mean- field result 
and the fiuctuations-corrected one are milder, since the 
gap is bigger. Once again, the gap is reduced by the 
fluctuations. Finally, figure © gives the specific heat 
through a larger temperature domain for the XXZ-like 
case. There is a rather large difference between the two 
theoretical curves (from mean-field and fluctuations) over 
the whole temperature range plotted. Comparing with 
the finite-size exact diagonalization (ECD)^^ results pre- 
sented, it is clear that the mean-field result has been 
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Specific heat vs. temperature tor {0.25, 1 , 1 ) at H 
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FIG. 4: Low-temperature specific heat for the XXZ case 
with anisotropy of 0.25, for different values of the external 
field. Four plots are given, corresponding respectively to H^ 
= 1.0, 1.4, 1.5 and 1.8. In each graph, the mean-field and 
fluctuations-corrected values are given by the two separate 
curves. The first three fields are below the critical field, the 
fourth somewhat above. Incorporating the fluctuations has 
reduced the gap. 



Specific heat vs. temperature for (1 , 0.33, 0.66): MF and MF+fluc 
H = 0.8 H = 1.2 




0.02 
T 



FIG. 5: Speciflc heat for the XYZ case, for two values of the 
magnetic field (below the transition). Both the mean-field 
result and the fluctuations-corrected one are plotted for small 
temperatures. 



greatly improved by incorporating fluctuations, and that 
at low temperatures {T/J < 0.2) the numerical curves 
tend to converge towards the theoretical ones when go- 
ing to higher size (for higher temperatures, the finite-size 
results are essentially those of the continuum limit) . 



C. Zero temperature excitation gap 

The mean-field prediction for the gap is easily found 
from the dispersion relation H1U|) . As mentioned previ- 
ously, below the transition two-particle excitations yield- 
ing incoherent spin-spin correlation functions dominate, 




FIG. 6: Specific heat for the XXZ-like case, plotted over a 
larger temperature range. Besides the mean-field and fluc- 
tuations results, exact diagonalization ones are presented for 
lattice sizes 14 and 16. 



and the gap is twice the single particle gap, whereas 
above the transition the one-particle excitation sets the 
gap. This yields 



OJo 



2IAU/1- 



/l2 



2i\t\ 
\h\ 



h\), 



\h\<\t~Ayti 

|t-AVt| < \h\ < \t\ 

\h\ > \t\ 



(40) 



To compute the gap in the presence of fluctuations, we 
numerically search for poles of the full Green's functions 
in the complex plane of ico. The Matsubara formalism ex- 
pressions we have obtained for general inverse tempera- 
ture p can be analytically continued away from the imag- 
inary axis of iuj when we formally take the limit of zero 
temperature, /3 ^ oo. We have developed a suite of nu- 
merical procedures to obtain the Green's functions in this 
limit. Poles of the full Green's function are then found by 
requiring explicitly that the denominator Ij^til) vanishes. 
The pole with the lowest real value of ioj then determines 
the position of the one-particle excitation gap. Explic- 
itly, we have written a numerical algorithm finding roots 
of 

a(fc, LUr)a{k, —LOr) — b(k, UJr)b{k, —LUr) = 0, (41) 

where 

a{ujr) = T.N{-(^r) + Go{k,U!r)/'DetQoik,UJr), 

b{uJr) = SsK) + ^o(fc, cj^)/Detao(fc, ^r)- (42) 

An effective dispersion relation ujr{k) is thus obtained, 
whose minimum can be found by numerically scanning 
along momentum values. The gap is then given by (Min^ 
here meaning the minimization of the function with re- 
spect to k) 



■'gap 



Miufc 2wr, \h\ < \t\, 
MiUfc LJr, \h\ > |t|. 



(43) 



Upon inspection of the results of these computations, 
some pleasant features emerge. The pure mean- field the- 
ory predicted a purely convex gap as a function of field 
below the transition; in the presence of fluctuations, this 
prediction gets greatly modified for fields below the crit- 
ical field. For example, in the case of figure ((HJ, the 
gap now reaches a maximal value somewhat below the 
transition, and decreases for decreasing field, therefore 
reproducing basic features expected of an exact solution 
of the problem. The zero-field gap (discussed in the next 
section) also agrees very well with the exact value: in fig- 
ure {Tj), the mean-field value at zero field was 22 percent 
away from the exact one, whereas fluctuations bring this 
discrepancy right down by two orders of magnitude to 
0.2 percent. In figure JHJ, the same comparisons yield 98 
percent versus II percent, again an increase in precision 
by about an order of magnitude. 



Gap vs. field tor (1 , 0.66, 0.33): DMRG, MF. and MF+tluc 
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FIG. 7; Gap as a function of field, computed in mean-field 
and with fiuctuations. The dot at zero field represents the 
Bethe Ansatz exact value. 



Gap vs. field for (1 , 0.33, 0.66): DMRG, MF, and MF+tiuc 



D. A note on the numerics 

We use the standard DMRG method with periodic 
boundary conditions'^ as an independent numerical tool 
to study the magnetizations and excitation gaps of the 
Hamiltonian Eq(I). We calculate the energies of the 
three lowest states of the systems, for system sizes up 
to N=100 spins. As expected, for H < He the gap 
between the lowest and the first excited state vanishes 
rapidly with growing system size reflecting the fact that 
the ground state is twofold degenerate in the thermody- 
namic limit. In Figs {THHJ we do not show these vanishing 
gaps but only display the energy diffrence between the 
second excited state and the lowest state. For a better 
comparison with the field theoretical results we extrapo- 
late our results using systems of size N = 30... 100 spins. 
As can be seen from Fig© our data for the gap extrap- 
olate well to the Bethe value, though the effects of the 
finite system size is still substantial even for 100 spins. 
The extrapolation procedure for non vanishing magnetic 
field H < He comes with an uncertainty which we es- 
timate to be 5 percent, this error is mostly due to the 
well known problem of unsystematic finite size effects in 
nonzero magnetic-fleld. 

For H > He the lowest gap is the energy difference 
between the lowest and the first excited state. Finite 
size effects in this region are negligible. 



IV. BETHE ANSATZ FOR ZERO FIELD 

The Hamiltonian of the zcro-flcld, fully anisotropic 
XYZ Heisenberg spin chain is related via trace identities 
to the transfer matrix of the eight-vertex model, solved 
exactly by Baxter—. Let us recall some aspects of the 
exact solution of the model, to which our perturbative 
results at zero field can be compared. For the zero-field 
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FIG. 8; Gap as a function of field, computed in mean-field 
and with fiuctuations. The dot at zero field represents the 
Bethe Ansatz exact value. 



anisotropic chain 

N 

^=E [J.s^s^^,+jysys]^,+j.s;s;^^ , (44) 

we consider the case Jz > Jj, > |Jic| > 0, which covers 
all parameter space for zero field. It is straightforward 
to map this restricted case to any of the cases we have 
considered above. 

Excitations over the ground-state are formed by 
spinous whose single-particle energy gap' is given by 



i^sp — t^,; 



sni2CK{u')u' 
2C ■ 



(45) 



Our notations are such that sn is Jacobi's sinus ampli- 
tudinus function of modulus k, where 



k = 



I 



l + V 






(46) 
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FIG. 9: Example of numerical data for the gap used in the 
extrapolation to A'^ — > oo, here for the zero- field XYZ magnet 
with exchange parameters (1,0.66,0.33). The circles repre- 
sent DMRG data for finite A'^, while the square is the exact 
zero- field value from the Bethe Ansatz. 



and sn/ is the function of modulus I. The parameter C is 
given by 



C = K{1) + 



2iri 

T+7' 



(47) 



where K{1) is the complete elliptic integral of modulus I, 
and T] is determined from the anisotropy parameters of 
the model according to 



sn^2r; = - 



V'-'z Jx ^ \ ' z '^y 



Jy Jx 



(48) 



Finally, the argument u and its complementary u' 
Vl — u^ are determined from 



c 



K{u) K{1')' 



(49) 



Since we are by definition below the critical field, 
spinous are always created in pairs, and therefore the 
physical excitation gap at zero field is given by 2Asp. 
We can plot the exact value for the gap from the 
Bethe Ansatz against the mean-field prediction, and 
the fluctuations-corrected prediction. We present three 
graphs of the gap as a function of the exchange parameter 
Jz, for three different choices of J^ and Jy. In figure (|10|l . 
we lie rather deep in the anisotropic limit, where the zero- 
field gap is rather large throughout the parameter range. 
The agreement between mean-field and the exact result is 
acceptable (the variation in the ferromagnetic end (neg- 
ative Jz) being 4.9 percent and that at the antiferromag- 
netic end being 8.6 percent), but incorporating fluctua- 
tions makes the agreement dramatically better (down to 
0.2 percent and 0.3 percent, respectively, so comfortably 



over an order of magnitude increase in precision) . In flg- 
ure l(TT)l . the x, y anisotropy is reduced, and therefore also 
the gap. This leads to greater disagreement, although the 
exact solution and fluctuations-corrected results are still 
extremely close. Here, the mean-field strays away from 
the exact solution by 15 percent and 43 percent for the 
ferromagnetic and antiferromagnetic ends, respectively, 
whereas the fluctuations results differ from the exact so- 
lution by only 1.8 percent and 2.2 percent respectively. 
Finally, figure ((T^ put us closer to the (gapless) isotropic 
case, and as expected shows the greatest deviations. For 
the purely antiferromagnetic case Jz > 0, however, the 
fiuctuations-corrected result reproduces the exact solu- 
tion to great precision up to Jz ~ 0.6. 

The Bethe Ansatz allows also the computation of ther- 
modynamic quantities in the model. The free energy at 
temperature T can be written 



/ = -TlnAo, 



(50) 



where Aq is the largest eigenvalue of the transfer matrix 
of the model. Its value is determined as 



Ao = 2[] 



Jz 



(1+P?) (51) 



where the pj parameters satisfy the coupled equations 

Jz 



2T 



G{pi) - Y. [<P{puPj) + 0(Pi,-p,)] = 27r(/ - 1/2), 



j=i 



G{p)=p 



^{p^ + Dl){p^ + Dl) 



(p2 + 1) 



0(p, ,) = 2 arctan I P^.A, + g(l +p^) ^^^^ 



^^{p^+Dl){p^+Dl)^ 



in which the anisotropy parameters appear as A^^ = 

Jx /Jz, A ^ = Jy/Jz, and D^ = \/l - A^, Dy = 

1 — A^. These equations can be solved numerically, 

which permits the computation of the temperature- de- 
pendent specific heat of the model at zero field. A com- 
parison of our perturbative results with the exact TBA 
result is given in figure ifT^ . 



V. CONCLUSION 

While the perturbative scheme that we have developed 
and implemented here does not constitute a definitive 
treatment of anisotropic Heisenberg chains in a field, 
which can only be achieved through a nonperturbative 
computation of Bethe Ansatz-type, it does provide an 
efficient and credible approach to the computation of 
approximate numerical values for physical quantities in 
general. The comparisons with DMRG and (for zero 
field) with the Bethe Ansatz that we have provided 



10 





Excitation gap at zero field for Jx = 


1,Jy = 0.3 




1.1 










1 


\>--, 

\>-, 




exact 

mean field 

fluctuations 




0.9 






- 


0.8 


\x_ 


- 


0.7 


^\^^ 


- 


0.6 




\: 


0.5 




'^-- 



Excitation gap at zero field for Jx = 1 , Jy = 0.9 



1 r 

0.9 - \ 

0.8 - 

0.7 - 

0.6 -"\ 

m 0.5 - 
0.4 

0.3 - 

0.2 - 
0, 




exact - 
mean field 
fluctuations 




FIG. 10: Excitation gap as a function of the exchange param- 
eter Jz, for strong x, y anisotropy. The straight mean- field re- 
sult lies up to 5 percent away from the exact solution, whereas 
the fluctuations-corrected result precisely overlaps it. 



Excitation gap at zero field for Jx = 1 , Jy = 0.5 
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FIG. 11; Excitation gap as a function of the exchange pa- 
rameter Jz, for weaker x,y anisotropy. Once again, the 
fluctuations-corrected result greatly improves the mean-field 
one. 



are a good guide to establishing the relative numerical 
precision to be expected from our theory, which ranges 
from very good for the "worst-case" scenario of a gapless 
XXZ chain, to essentially indistinguishable from the 
exact solution (errors of order 0.2 percent) for spin 
chains with strong structural XYZ anisotropy. In broad 
terms, an increase in precision of at least an order of 
magnitude is achieved by incorporating the fluctuations 
as we have done here, which makes our approach good 
enough for many practical purposes. Here, we have 
concentrated on equilibrium thermodynamic quantities 



FIG. 12; Excitation gap as a function of the exchange param- 
eter Jz, for weak x,y anisotropy. The fluctuations-corrected 
result greatly improves the mean-fleld one, but the agreement 
is less spectacular than in the above two cases, due to the rel- 
ative smallness of the zero-field gap for weak anisotropy. 
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FIG. 13; Specific heat for an XYZ magnet in zero field, com- 
puted from the Bethe Ansatz and compared to the straight 
mean-field and fluctuations results obtained in two different 
ways (interchanging the Y and Z anisotropies, i.e. our per- 
turbation parameter). The fluctuations result for Jz = 0.33 is 
indistinguishable from the exact result on the resulting plot. 



like magnetization, and basic dynamical-related ones 
like the excitation gap. It would be of great interest to 
extend our treatment to purely dynamical ones like the 
spin structure factor, which we will do in a forthcoming 
publication. 
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Specific heat vs. temperature at H = 0: Exact Complete Diagonalisation and TBA 
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